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Abstract 

We construct new electrically charged solutions of the Einstein-Maxwell equations with negative cos- 
mological constant in general odd dimensions d = 2n + 3 > 5. They correspond to higher dimensional 
generalizations of the squashed black hole solutions, which were known so far only in five dimensions. 
We explore numerically the general properties of such classical solutions and, using a counterterm pre- 
scription, we compute their conserved charges and discuss their thermodynamics. 

1 Introduction 

Motivated by recent advances in String/M- Theory (as well as theories with large extra dimensions in TeV 
gravity), black holes in higher dimensions have been actively studied in recent years with many surprising 
results (for reviews of this topic see e.g. [1, 2]). In d = 5 spacetime dimensions, an interesting class of 
solutions is provided by the so-called squashed Kaluza-Klein (KK) black holes, whose horizon geometry 
is a squashed three-sphere [3, 4, 5, 6]. Their geometry is asymptotic to a non-trivial S 1 bundle over a 
four-dimensional asymptotically flat spacetime, which is also the asymptotic geometry of the Kaluza-Klein 
magnetic monopole [7, 8] . Such black holes look five-dimensional in the near- horizon region, while at infinity 
(asymptotically) they look like four-dimensional objects with a compactificd fifth dimension. KK black hole 
solutions in the presence of matter fields, are generally found by solving the Einstein equations by brute 
force. For instance, a solution describing a static KK black hole with a U(l) electric field has been found 
in [9], the corresponding non-Abelian generalization being described in [10]. Remarkably, with hindsight, 
many such KK solutions can be generated from known solutions by applying a 'squashing' transformation 
on suitable geometries [11, 12, 13, 14, 15, 16]. 

Given the existence of such exotic objects in asymptotically flat settings, it is naturally to ask if such 
solutions do have counterparts in presence of a cosmological constant. While on physical grounds one would 
expect such objects to exist, so far exact solutions have proven to be elusive and one is once again compelled 
to recourse to numerical methods for solving Einstein's field equations. Indeed, finding higher dimensional 
solutions in (anti)- dc Sitter ((A)dS) back backgrounds is even harder than the asymptotically flat case, as 
most of the solution-generating techniques break down in presence of a cosmological constant. Nonetheless, 
there are some approximate solutions and also numerical ones that show the existence of squashed AdS black 
holes in five dimensions [17, 18]. Futhermore, extremal charged squashed black holes have been studied in 
presence of a positive cosmological constant in [19] . 



The main purpose of this work will be to extend sonic to these results to higher than five dimensions 
and to construct new solutions of the Einstein-Maxwell theory that describe charged squashed AdS black 
holes in general odd dimensions. Our construction is based on the simple observation that one can write 
the metric of an odd-dimensional (round) sphere S 2n+1 as an S 1 fibration over the complex projective space 
CP n : 

cM 2 2n+1 = (dz + Anf + dSl (1) 

where cE 2 is the metric on the 'unit' CP n space, A n its Kahler form and <Kl\ n+l is the metric on the 'unit' 
round sphere S 2n+1 (we are using the results and conventions from [20]). Using this ansatz, the squashing of 
the sphere S 2n+l will be described by a squashing parameter N that will multiply the Kahler form A n in (1). 
Then, in the un-squashed case (here N = 1), our solutions would correspond to the usual Schwarzschild-AdS 
black holes with spherical horizons, while for N = they correspond to a version of the odd-dimensional AdS 
black strings found in [21, 22, 23]. For general values of N we shall show that there exists a class of black holes 
in AdS-backgrounds whose horizons correspond to a squashed S 2n+1 sphere. Since we were unable to find 
exact solutions describing such squashed AdS black holes, the configurations in this paper are constructed 
numerically by matching the near-horizon expansion of the metric to the asymptotic Fcfferman- Graham 
form. 



2 The model 

2.1 The ansatz and the equations of motion 

We start with the Einstein-Maxwell action principle in d-spacetime dimensions: 

d d xy/^g(R — 2A - F 2 ) — / d d ~ x x^fK, (2) 



16nGd Jm &nGd Jqm 

where Gd is the gravitational constant in d dimensions, A = — (d— l)(d— 2)/(2£ 2 ) is the cosmological constant 
and F = dA is the electromagnetic field strength. Here Ai is a d-dimensional manifold with metric . K is 
the trace of the extrinsic curvature K a b = — 7^V c n(, of the boundary dM with unit normal n a and induced 
metric j ab . 

By setting the variations of the action (2) to zero, one obtains the Einstein-Maxwell system of field 
equations: 

G,u- {d ~ 1 ^~ 2) 9,u = 2T, v , V^=0, (3) 

together with a Bianchi identity for the electromagnetic field dF = 0. Here G^ v is the Einstein tensor and 
T M „ is the stress tensor of the electromagnetic field. 

We consider here the following parametrization of the d-dimensional line element (with odd values of 
d = 2?i + 3 > 5): 

ds 2 = a(r) (dz + NA n f + + ^—^-r 2 dY,l - b(r)dt 2 , (4) 

f{r) [d-4) 

where the (d— 3) -dimensional metric gE 2 is the metric on the 'unit' CP n space, A n is its Kahler form and 
N is the squashing parameter. The direction z is periodic with period L = 27rjv v / ( d ^ w here k is an 

integer number. Here the unusual normalization in front of the CP n metric dY? n has been chosen in order 
to make contact in the case N = to the results obtained previously for the charged black strings in AdS 
backgrounds [23], as the normalized metric ^~^ r 2 rfE% has the same normalization as the unit sphere S d ~ 3 . 
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We shall also consider an electric ansatz for the electromagnetic field, 

A = V{r)dt. (5) 

The Einstein equations with a negative cosmological constant imply then that the metric functions a(r), b(r) 
and f(r) are solutions of the following equations: 

, 4(d-4) 2 a 2 , Jo! , b'\ 2(d-l)r 2(d-4)(l-/) 4r/ . 



f + — '- — N z + f 1 — H - — V u = 

1 + (d-l) 2 r3 +/ \a b) P r + (d - 2)6 

, _ (rf-3)(rf-4)6 _ 1 (d-l)(d-4)6 _ (d - 3)6a' (rf-4)6'(l-/) (d-l)^ 

r 2 1 / j Pf ra rf Pf 

a'b> b' 2 (^4) 2 a((d-3)6 + 2r6Q A 2 /3d -8 rb> \ 2 _ 

2a 6 (d-l) 2 r 4 / V d - 2 + (d - 2)6/ ' W 

1 / 2a6 



((d - 3)(d - 4)£ 2 {f - 1) - (d - 2)(d - l)r 2 ) + 2(d - 3)a6' 

N 2 +4raV' 2 ) = 0. 



rb' + 2{d-3)b\rf£' 

2(d-4) 2 (d-3)a 2 6 



(d-l) 2 r 3 / 

One can then easily solve the Maxwell equations in (3) by taking 



,d-4,i (d _ 3 ) f r 1 / 6(0 



where Q, $ are integration constants, which will be later related to the electric charge and chemical potential, 
respectively. After replacing the corresponding expression of V in (6), the problem reduced to solving a 
system of three coupled ordinary differential equations, with given parameters N, I and Q. 

2.2 Asymptotics 

The corresponding asymptotic expansion of the metric functions for a number d = 2n + 3 of spacetime 
dimension is given by 1 

«« - ^ + ( E /2 «,(;) 2j +c a iog(^)(^) d - 3 +C2 (^r 3 + o ( ^), 

3=0 
2 (d-B)/2 

*(r) - ^+ E ^^) 2j +C h log(J)(^) d - 3 + c t (^)^ + 0(^), (8) 
j=o 

2 (d-5)/2 

/m = r -,+ E /,A 2j +2c / io4)A- 3 + (c , +Ct+C0 )A- 3 + o(^)- 

The coefficients a,, 6j, /, as well as £ a , C/ an d c o depend only on d, iV, £ and have a complicated expression, 
with no apparent pattern, except for the lowest order, where one finds: 

d-4/ d-4N 2 \ d-4/ (d-4) 2 iV 2 \ a 

fl0 = d^V-d^W)' b °=d^{ 1 -Td--W2p)> ^ /0 = &0+ d^2- (9) 



iNote that in all (odd) spacetime dimensions, one finds the asymptotic expression of the Riemann tensor R^" = — (<5^<5JJ 
<5^)/£ 2 + 0(l/r d -4). 
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To give a flavour of these expressions, we give here the asymptotic form the solution for d = 7 only 
(similar expressions were found for d — 9, 11, while the corresponding relations for d = 5 were given in [18]) 

t , . r 2 9 iV 2 U 2 / 17iV 2 \ / N 2 \ £ 4 / 9 
f{r)= P + W 5^ + 320^ l^^J ( V 2 -^J + ^( V 1600 +Ct + C2 

9 ?u2 ( og i 173jV 2 94W 4 ^ , . „ N , 

3iV [Sb+—p 9£ 4 / N 2 \ ( 25N 2 34iV 4 \ r „Aogr, 



6400£ 2 J 1600r 4 V" P J \® + I 2 £ 4 ) 1 ° g( £ ) + ° ( r 6 ^ 

, . 7* 2 3 3iV 2 U 2 / 43iV 2 \ / 7V 2 \ ^ 4 

aW = ^ + 4-^ + 64(p( 6 -^ 2 -j( 2 -^J + ^ (10) 

3<? 4 / N 2 \ ( n 97iV 2 1537V 4 \ ,7\ ^,logr. 
'2--^ 18+-= -n- logb) + 0(-^), 



3200r 4 V ^ 2 / V « 2 ^ 4 / « r 

, . r 2 3 N 2 U 2 / 1EV 2 \ / N 2 \ t 

h ^ = p + l-W + M^\ ~^) \ 2 -l 2 -) + ^ Ct 

U A / N 2 \ ( 53iV 2 517V 4 V ,r, 4og 



18+^5 33- lo g(7) + 0(-S-) 



[ r 



3200r 4 V I 2 ) \ i 2 £ A J oy £' v r 6 

For any (odd) value of d, terms of higher order in l/r depend on the two constants ct and c z and also on 
the charge parameter Q (which enters only the sublcading order terms in (8)). The asymptotic expression 
of the electric potential is V{r) = $ - Q(jEj )^ (d ~ 3) {d J 3 f rd - 3 +■■■■ 

The constants (ct, c z ) are found numerically starting from the following expansion of the solutions near 
the event horizon (taken at constant r = rjj) and integrating the EM equations towards infinity (this 
expansion corresponds to a nonextremal black hole): 

a(r) = a h + oi(r - r H ) + 0{r - r H ) 2 , b{r) =h{r- r H ) + 0{r - r H ) 2 , f{r) = fi{r - r H ) + 0(r - r H ) 2 . (11) 

The only free parameters in the event horizon expansion are ah, b±. One finds e.g. 

f- 1 (i(d 4) ° 2 | W-^h 4 a ,(^4) 2 iV 2 \ 

h ^ { 2(d 4) " a h (d-2) r 2 ^ + — P (d^T^ ) ' (12) 

n ( d— 4 \et— 3 



2a h d-2 

Ol 



Q 2 , ahjd-l)^ al(d-4Y(d-3) iV 2 

2(d-4) + p + {d-^ 1 ^ 



r H (A _ A\n _ 2 (^t)"' 3 Q 2 , ^(d-l)r% _ 2a 2 (d-4)2(rf-3) AT2 ' 

the higher order therms being very complicated, and, for simplicity, we shall not list them here. 

Let us also note that in the near horizon region one finds the following expansion of the electromagnetic 
potential, 

V(r) — Vq + fgEi) 2 '\f^Jl^ r ~ t h) + ~ Th ) 2 ' wnere is an integration constant. Notice 

now that one can always set Vq = such that V(rjj ) = 0. Then the physical significance of the quantity $ 
in (7) is then that it plays the role of the electrostatic potential difference between the infinity and horizon. 

The condition for a regular (and nonextremal) event horizon is f'{rh) > 0, with b'(rh) > 0, which, for 
given th, imposes an upper bound on Q, 

^ < - 2 , (^f V* (<- « + « - # - -»^f f ) ■ m 

2.3 Global charges and thermodynamical quantities 

The computation of the mass/energy of these solutions is basically similar to that performed in [22] for the 
black string limit of these solutions. Employing a a counterterm prescription, the total action is supplemented 
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with two boundary terms, I = Iq + 7° t + I* t , with the standard AdS counterterm action [24, 25]: 

d-i„r^f d - 2 ^- 4 ) D e 3 e(d~6) / DQb d-i d2 



8^G d 7 V 7 1 £ 2(d-3) 2(d-3) 2 (d-5) ^ ab 4(d-2) 

^ 5 e(d-8) ^3d-i 6 d 2 -i D3 

-KK K a fc — . — r-^ K 



(d - 3) 3 (d - 5)(d - 7) \A(d - 2) aD 16(d - 2) 2 
-2R ab R cd R acM -_l_Lv a RV a R + V c R ab V c R a ^ + ...| , (14) 

where R and R ah are the curvature and the Ricci tensor associated with the induced metric 7. The series 
truncates for any fixed dimension, with new terms entering at every new even value of d, as denoted by the 
step-function (0 (x) = 1 provided x > 0, and vanishes otherwise). 

The extratcrm 7^ t is due to the presence of \og(r/£) terms in the asymptotic expansions of the metric 
functions (with r the radial coordinate) [26]: 



128 



RR ab R Qb - ^R 3 - 2R ab R cd R acbd -^R a6 V a V t R + R a6 DR afc - S d , 7 + . . . J 



Using these counterterms in the total action, one computes the boundary stress-tensor: 

2 SI 



ab 



/— 7 £7 

Following the usual approach, we consider now a standard ADM decomposition of the metric on the boundary, 
- 1ab dx a dx b = ~V 2 dt 2 + a ij (dy i + V l dt)(dy j + V j dt), 

where V and V 1 are the lapse function, respectively the shift vector, and y l , i = 1, . . . , d — 2 are the 
intrinsic coordinates on a closed surface £ of constant time t on the boundary. 

Then a conserved charge 0^ = ^ s d d ~ 2 yy / au a £ i b T a i ) , can be associated with the closed surface £ (with 
normal u a ), provided the boundary geometry has an isometry generated by a Killing vector £ a . 

For the line element (4), the conserved mass/energy M is the charge associated with the time trans- 
lation symmetry, with £ = d/dt. There is also a second charge associated with the compact z direction, 
corresponding to the squashed black hole's gravitational tension T . 

The computation of the boundary stress-tensor T ao is straightforward and we find the following expres- 
sions for mass and tension: 

M = Mo + M<?> , with M = i-1 [c, - (d - 2)ct] LV n , (15) 

167rG 

T = To + T c {d) , with To = -— [(d - 2)c 2 - ct] V„ , (16) 

167rG 

where V n is the total area of the 'unit' CP n . Here Af c (d) and 7; (d) are Casimir-likc terms, 



M ( d ) * 4 v n L 



T ( d) = 



16ttG 



16ttG 



1 TV 2 o 1 A 2 A 4 A 6 

12 (1 - 16?^ " ^OO^ 1008 + 40 V ~ 26 V + + ■ • 

1 19A 2 (32/ 2 - 7N 2 ) 1 TV 2 A 4 A 6 



The Hawking temperature as obtained from evaluating the surface gravity or requiring regularity on the 
Euclidean section is: 

1 



Th = ^VKfi, (17) 
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with fi given by (12). The area Ajj of the black string horizon is given by 



Ah = r*- d V d - 3 L^ . (18) 

As usual, one identifies the entropy of the charged squashed black hole solutions with one quarter of the 
event horizon area 2 , S = An/^Gd- 

Finally, for a charged solution, the electric field with respect to a constant r hypcrsurface is given by 
E p = g^ p F p „n". The electric charge Q of the charged solutions is computed using Gauss' law by evaluating 
the flux of the electric field at infinity: 

2 = ^//-W^ = -^. (19) 

If Aft is the electromagnetic potential, then the electric potential $, measured at infinity with respect to the 
horizon is defined as: 

$ = ^1^00-^x^1^, (20) 
with a Killing vector orthogonal to and null on the horizon. 



3 The solutions 

Solutions of the system (6), (7) with the asymptotic behaviour (8), (11) are known to exist in two different 
limits. First, for a vanishing squashing parameters, N = 0, one recovers the AdS black strings 3 discussed 
in [21], [22]. Their generalization with an electric U(l) field has been considered in [23]. The second known 
configuration is found for N = No = \J jE^^ an d corresponds to the RNAdS black hole (with a nonstandard 
parametrization of the spherical sector). For our purposes it is helpful to look at this solution from the 
viewpoint of the shooting procedure. Assuming that there is an horizon at r = th , its expression within the 
parametrization (4) has a(r) = and 

m . m _ £ (1 _ cf?-i) + (i - e>«) + <i§§^ (i - ^ p., 

the global AdS spacetime being recovered for Q = rn = 0. 

An analytic solution of the equations (6) for N ^ No appears to be intractable, and one is contrained to 
consider numerical methods. In our approach, the equations (6), (7) are treated as a boundary value problem 
for r € [rn , oo] (thus wc did not consider the behaviour of the solutions inside the horizon) . There we employ 
a collocation method for boundary-value ordinary differential equations, equipped with an adaptive mesh 
selection procedure [27]. The mesh include 10 3 — 10 4 points, the typical relative accuracy, as estimated by 
the solver being of order 10 -8 . The problem has four input parameter: r#, N, £ and Q. The event horizon 
data a(rn), b' '(rsf) and the coefficients at infinity q, c z and $ are read from the numerical output. 

We have found in this way charged squashed black hole solutions with AdS asymptotics for d = 5, 7 and 
9, the first two cases being studied in a systematic way. Thus they are likely to exist for any d = In + 3 > 5. 
For all the solutions we studied, the metric functions a(r), b(r), f(r) and the electric gauge potential V(r) 
interpolate monotonically between the corresponding values at r = and the asymptotic values at infinity, 
without presenting any local extrema. 



2 This relation has been derived in [22] by using Euclidean techniques. This is not possible for the electrically charged 
solutions in this work, which do not solve the equations of motion for an Euclidean signature. However, it can be deduced using 
the quasi-Euclidean approach. 

3 Note that these solutions are not known in closed form; however, they were recovered independently by several different 
groups, with different numerical methods. 
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Figure 1. The profile of a typical charged, squashed d — 7 black hole (left) is shown together with a that of a 
d = 7 soliton (right). 

As a typical example, in Figure 1 (left) the metric functions a(r), b(r) and f(r) as well as the electric 
gauge potential V(r) are shown as functions of the radial coordinate r for a d = 7, N = 0.7 solution with 
rjj = 0.5, Q = 0.3, I = 1. One can sec that the term r 2 /£ 2 starts dominating the profile of a(r), b(r) and 
f(r) very rapidly, which implies a small difference between the metric functions for large enough r, while 
the gauge potential V(r) approaches very fast a constant value. 

In our numerics, we have fixed the length scale by setting i = 1 and study the dependence of the solutions 
on th, N and Q. Also, since only N 2 appear in the equations, wc shall restrict to positive values of N. 

Starting with the dependence of the solutions on the event horizon radius rn for fixed squashing and 
electric charge, our numerical results show that the picture valid for the RN-AdS solution is generic. For 
example, for large enough Q, one finds only one branch of thermally stable black holes, i.e. their entropy 
increases with the temperature. Also, for Q ^ 0, the solutions exist for a ru > ro, with ro a solution of the 
equation 



2 N 2 _ (d-l) d - 1 (d-2) /ro x 2(^-3) / r (d_ 3 )(d-4) 

2(d-4) d - 2 \£J [} l> {d-1) 2 



= 0. (22) 



As Th — > ro, an extremal black hole with Tjj = is approached, the corresponding near horizon solution 
being 



where 



(d-3)(d-4) 



(^2)(^l)(# + ^f 



V(d-2)(d-l)J^+ ^^" (24) 



In Figure 2, we show the dependence of the event horizon area, mass and tension on the Hawking temperature, 
for three different values of Q. These plots retain the basic features of the solutions we found for other values 
of d, N (note that in this paper we set L = V n /16irGd = 1 in the numerical values for the mass, tension and 
charge; also we subtracted the Casimir energy and tension). 
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Figure 2. The mass M, tension T and event horizon area Ah are shown as function of the Hawking temper- 
ature Th for d — 5 squashed black hole solutions with three different values of the electric charge parameter Q. 
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Figure 3. A number of parameters are shown as a function of the electric charge (left) and of the squashing 
parameter N, for d = 7 squashed black holes with fixed event horizon radius. 



The case Q = is special, since the solutions exist for any rjj > 0. As expected, the properties in this case 
are rather similar to those of the Schwarzschild-AdS black hole. For example, for any N, their temperature 
is bounded from below, and one finds two branches consisting of smaller (unstable) and large (stable) black 
holes. 

The Q = solutions with a vanishing event horizon radius, r# = 0, are of special interest, corresponding 
to a deformation of the global AdS background. These are soliton with a regular origin, the behaviour of 
the metric functions as r — > being 

«M = |^|^ + 0(r 4 ), Kr)=bo+ b -f^f + 0(r% f(r) = ^ + f 2 r> + 0(r% (25) 

in terms of two parameters 6o, For large r, they also approach the asymptotic form (8). 

Interestingly, an approximate analytic expression can be constructed in this case, viewing the solutions 
as perturbations around the global AdS background. The solutions are constructed order by order, the 
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perturbation parameter being e = N/N — 1. However, the resulting expressions are exceedingly complicated 
(and not very enlightening) except for the hrst order in e. Then one finds the following distorsion of the 
global AdS background 

r 2 1 / 5 I 2 ( 1 i 2 4r 2 \\ 

/(') = p + J + « (-g + si (-1 + 4( 5 + j^Hogd + , (26) 

2 ^ 2 . , 4r 2 A ., , r 2 ( 1 € 2 , , 4r 2 



for d = 5, 



fl M = 72 + M - 1 + 4^2 lo sd + i2-))> 6 M = 72 + e ( - 2 + 12^2 + it) 



„ . r 2 1 / 4 2« 2 / 5£ 2 6^ 2 5£ 2 . , ,„ 2r\\\ 

^ r2 f 3 3£ 2 / £ 2 , „ 2r\\\ t/ . r 2 1 ( \ U 2 / ^ 2 , ^ 2r 2 . 
a « = F + %-2 + 2^( 1 -2^ l0g(1 + ^ ) JJ' 6(r) = F + 2 +e (-2 + l^ 

for d = 7, its expression for higher d becoming too complicated to display it here. The nonperturbative 
solutions with arbitrary N are constructed again numerically, a typical soliton profile being exhibited in 
Figure 1 (right). 

Returning to solutions with an event horizon, we show in Figure 3 the dependence of the d = 7 squashed 
black holes on the electric charge (left) and on the squashing parameter N (right), for a fixed event horizon 
radius. One can see for example that, for given r#, the mass, tension and entropy of solutions decrease 
monotonically with N. Nontrivial solutions appear to exist for all values of N, the value A*o = 
separating prolate metrics from the oblate case (N > No). 



4 Further remarks 

The specific class of solutions that we constructed here are likely to be relevant in connection to the so- 
called 'fragility' of the AdS black holes. As described in [28] (see also [29]) it turns out that AdS black 
holes can become unstable to stringy effects when their horizon geometries are sufficiently distorted. By the 
AdS/CFT correspondence, this implies that their field theory on the boundary can also become unstable. In 
five dimensions this phenomenon was explicitly checked in [17, 18]. More specifically, if the spherical horizon 
of a black hole is sufficiently squashed then the theory on the boundary ceases to be stable since the scalar 
sector of the theory becomes tachyonic for sufficiently large squashing. At first sight, the squashed AdS 
black hole geometry does not exhibit any signs of gravitational and/or thermodynamic pathologies for this 
same amount of squashing of the black hole horizon that causes the instability of the boundary field theory. 
However, as shown in [28], the source of the instability in the background of the AdS squashed black hole is 
related to the stringy pair-production of branes in this geometry, a phenomenon discovered by Seiberg and 
Witten [30]. 

The Q = configurations we described in this paper reduce in five dimensions to the squashed black 
holes studied in [17, 18], while providing the needed generalizations to higher than five odd dimensions. 
As such, they will be directed relevant in connection to the fragility of the squashed AdS black holes in 
higher than five dimensions. In particular, we expect that for certain values of the squashing parameter 
the squashed geometry becomes unstable to the phenomenon of pair-production of branes and this will 
signal potential instabilities of the field theory on the boundary. Let us consider as an example the seven- 
dimensional squashed geometry. In this case the Ricci curvature scalar of the rescaled boundary geometry 
has the expression: 

ft = W 
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It is now clear that for sufficiently large values of N the Ricci scalar curvature can be constant and negative. 
This is a sufficient condition for the onset of the non-perturbative Seiberg-Witten stringy effect of pair- 
production of branes in this geometry. According to the AdS/CFT conjecture, this signals an instability in 
the boundary conformal field theory. 

As avenues for further research, let us mention that we have also found solutions that correspond to 
black hole geometries with squashed horizons in absence of the cosmological constant. These solutions have 
a different asymptotics expansion than (8) and will be discussed elsewhere. Also, one should be able to 
further generalize these solutions to include the effects of rotation. However, as these new solutions deserve 
separate studies, we shall report them in further work. 
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